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We present a quenched lattice calculation of all six form factors: vector [fi(q 2 )], weak mag- 
netism [f2(q 2 )], induced scalar [f3(q 2 )], axial-vector [<?i(? 2 )], weak electricity \g2(q 2 )], and induce 
pseudoscalar [<73(5 2 )] form factors, in hyperon semileptonic decay H° — » E + Zz7 using domain wall 
fermions. The q 2 dependences of all form factors in the relatively low q 2 region are examined in 
order to evaluate their values at zero momentum transfer. The H° — > E + transition is highly sen- 
sitive to flavor SU (3) breaking since this decay corresponds to the direct analogue of neutron beta 
decay under the exchange of the down quark with the strange quark. The pattern of flavor SU (3) 
breaking effects in the hyperon beta decay is easily exposed in a comparison to results for neutron 
beta decay. We measure Sf/(3)-breaking corrections to /i(0), /2(0)//i(0) and gi(0)//i(0). A sign 
of the leading order corrections, of which the size is less than a few %, on /i(0) is likely negative, 
while /2(0)//i(0) and <?i(0)//i(0) receive positive corrections of order 16% and 5% respectively. 
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The observed pattern of the deviation from the values in the exact 517(3) limit does not support 
some of model estimates. We show that there are nonzero second-class form factors in the H° — > E + 
decay, measuring /3(0)//i(0) = 0.14(9) and <?2(0)/<?i(0) = 0.68(18), which are comparable to the 
CNJ ■ size of first-order SU(3) breaking. It is also found that the SU(3) breaking effect on <?3(0)/(?i(0) 

agrees with the prediction of the generalized pion-pole dominance. 
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I. INTRODUCTION 

The latest lattice calculations of K — » irlu (Kis) semileptonic decays have been greatly developed with high 
precision [2,12)01. The results for the vector form factor /+(0) of Ki 3 decays, which deviates from unity due to flavor 
SU(3) breaking, can be used to provide a very precise determination of the element V us of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [1, Q. Those first-principles calculations contribute greatly to a stringent test of the CKM 
unitary through the first row relation |V u d| 2 + |Ki S | 2 + \V u b\ 2 = 1 [f|. On the other hand, AS = 1 semileptonic hyperon 
decays provide alternative determinations of \V US \ The consistency of the values of \V US \ determined from different 
experiments needs to be confirmed. As we will explain later, however, the determination of |V^ S | from the hyperon 
decays should be affected by larger theoretical uncertainties than those of Ka decays [§| . This is simply because the 
lack of a reliable theoretical calculation of the leading symmetry-breaking corrections in the hyperon beta decays. 
Indeed, the conventional Cabibbo model where flavor SU(3) breaking effects are ignored, is commonly used for 
the analysis of the hyperon beta-decay data [2|. 

An essential difference from the case of K13 decays is that the axial- vector current also contributes to the transition 
in the hyperon beta decays. Thus, the precise determination of \V US \ in the hyperon beta decays requires information 
of the ratio of the axial- vector to vector form factors gi(0)//i(0) in addition to the vector form factor /i(0) [j|. Here 
we recall that the various ratios g±(0)/fi(0) in the hyperon beta decays also provide vital information to analysis 
of strange quark spin fraction of the proton spin, together with the polarized deep inelastic scattering data [lfj. 
However, such analysis heavily relies on the Cabibbo model. The flavor SU(3) breaking introduces the systematic 
uncertainty on the strange quark contribution to the proton spin. This issue is still under debate [12] • Recently, a 
new analysis, where the first-order corrections of both SU (3) and SU (2) symmetry breaking are properly taken into 
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account within the Cabibbo model, was proposed by Yamanishi [12j. The author has reported that flavor breaking 
effects significantly affect the evaluation of the amount of spin carried by the strange quarks inside the proton [12| . 

In the hyperon decays, a model independent evaluation of St/(3)-breaking corrections is highly demanded for both 
the CKM unitary test and the proton spin problem. The weak matrix element of the hyperon beta decays can be 
calculated with high accuracy from first principles using the techniques of lattice QCD, similar to what is achieved in 
the case of decays [l|, [HQ. There is a single lattice study to be completed for a specific hyperon decay, E~ — > nlv. 
The lattice simulation have been performed by Guadagnoli et al. with 0(a)-improved Wilson fermions in the quenched 
approximation [l3[ [57j |. As we are interested primarily in flavor ST/ (3) breaking effects not in the absolute values of 
each form factor, we choose a different hyperon decay, S° — » in this study. This is simply because the S° — » E + 

transition is the direct analogue of neutron beta decay under the exchange of the down quark with the strange quark. 
The flavor SU (3)-breaking can be easily exposed through a comparison with results of neutron beta decay. In Ref . [l4| , 
for the axial-vector coupling of neutron beta decay, gi(0) = 1.212(27) in the chiral limit is obtained from quenched 
lattice QCD calculations with domain wall fermions (DWFs). It underestimates the experimental value of 1.2695(29) 
by less than 5%. Other relevant weak form factors have been subsequently investigated in Ref. [15]. We naturally 
extend the quenched DWF calculation for investigating SU (3)-breaking corrections to the hyperon semileptonic decay 
form factors. 

Our paper is organized as follows. In Sec. [TTJ we first present a brief introduction of the hyperon beta decays. In 
Sec. [m details of our Monte Carlo simulations and some basic results are given. We also describe the lattice method 
for calculating the baryon beta-decay form factors. Section IIVI is devoted to our determination of the scalar form 
factor fs(q 2 ), which will be defined in the next section, at finite momentum transfer. We describe the interpolation 
of the form factor to zero momentum transfer and the chiral extrapolation in order to evaluate the value of /i(0) at 
the physical point. We also present our estimate of \V US \. Section fVl gives the result of the ratio of ?i(0)//i(0) at the 
physical point. We then discuss flavor St/ (3)-breaking effects appeared in <7i(0)//i(0) for the S° — > E + transition in 
comparison to neutron beta decay. The results for the other form factors including the second-class form factors fz 
and (?2 are presented in Sec. IVIl Finally, in Sec. IVII1 we summarize our results and conclusions. 



II. HYPERON BETA DECAYS 



The general form of the baryon matrix element for semileptonic decays B — > blv is given by both the vector and 
axial- vector transitions: 

(b(p')\V a (x) + A a (x)\B(p)) = u b (p') {Ol{q) + Ot{q)) u B ( P y-*, (1) 

where q = p — p' is the momentum transfer between the initial state (B) and the final state (b) which belong to the 
lightest J p — 1/2+ SU(3) octet of baryons (p, n, A, S,S). The vector and axial- vector currents are defined as V a (x) = 
u(x)j a d(x) and A a (x) = u(x)^ a ^/^d(x) for AS = decays, and V a (x) = u(x)j a s(x) and A a (x) — u(x)j a j^s(x) for 
AS = 1 decays. Six form factors are needed to describe the hyperon beta decays: the vector (/i), weak magnetism 
(f2) and induced scalar (fa) form factors for the vector current, 

O a (?) - 7 «A (I ) + ^u^M h + ^J^TIh (2) 
and the axial-vector (gi), weak electricity (52) and induced pseudo-scalar (173) from factors for the axial current, 

O a (?) = 7«75 5l (? ) + ^/W^—^ + ^75^—^-, (3) 

which are here given in the Euclidean metric convention (we have defined a a p — ^[70,7/3]-) [5Sj |. Here, Mb (Mb) 
denotes the rest mass of the initial (final) state. Note that although the sign convention of the /a and (73 form factors 
is opposite in comparison with that of Ref. [l6j , in our convention both gi and 53 form factors are positively defined 
for neutron beta decay. In addition, our adopted normalization of X/(Mb + M b ) : instead of 1/Mb that adopted in 
experiments, is theoretically preferable for considering the time-reversal symmetry on the matrix elements [591 ] . 

For convenience in later discussion, we consider the scalar form factor fs(q 2 ) for the vector-current form factors 
given in Eq. ([2]): 

/#-V) - /f- V) + M^TM-Jf ^V), (4) 
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which can be defined through the matrix element of the divergence of the vector current as (b(p')\d a V a (0)\B (p)) = 
(Mb — MB)fs{l 2 )ub{p')uB{v) (see Appendix A for details), and also introduce a particular linear combination of the 
axial-vector-current form factors given in Eq. (|3|) as 



gf - V) = 9? ^ V) - 



M B - Mi 
M B + M b 



V), 



(5) 



which is defined in an alternative parametrization of O^(q) (see Appendix B for details). Both f$ and <?i form factors 
are relevant in lattice calculations (l3| . 

In the literature, the vector and the axial-vector form factors at zero momentum transfer are called the vector 
coupling gy = /i(0) and the axial- vector coupling gj± — <?i(0), respectively. According to Weinberg's classification [13], 
the terms /a and 172 are known as the second-class form factors, which are identically zero in the certain symmetric 
limit (iso-spin symmetry, [/-spin symmetry or U-spin symmetry as SU (2) subgroups of the flavor SU (3) symmetry) 
within the standard model. For an example, the second-class form factors in AS = decays such as neutron beta 
decay are prohibited from having nonzero values because of G-parity conservation in the iso-spin symmetry limit 
(m u = rrid) @- For AS = 1 decays, the U-spin symmetry (m,j = m s ) plays a similar role instead of the iso-spin 
symmetry. Observation of nonzero second-class form factors corresponds to the direct signal of flavor SU (3)-breaking 
effects in the hyperon beta decays. On the other hand, the terms f% and 53 are suppressed in the evaluation of the 
beta-decay transition amplitude by a factor [mi /(Mb + Mb)] 2 where mi is the charged lepton mass. Therefore, for 
the decay B — > bev e , their contributions can be safely ignored. As no accurate experiment has yet been performed on 
muonic hyperon decays, it is hard to access information of f 3 and 173 form factors in present experiments. 

In the exact SU(3) limit, the vector couplings are simply given by SU(3) Clebsch-Gordan coefficients as /i(0) = 
fkim , while the axial- vector couplings g\ (0) are governed by two parameters F and D as g\ 

(F = 0.475(4) and D = 0.793(5) are quoted in Ref. [l2| for the conventional Cabibbo fit). Here, dkim is the totally 
symmetric tensor of the SU(3) group. Moreover, the conserved vector current (CVC) hypothesis becomes valid in 
this limit. The vector part of the weak current is a conserved current like the electromagnetic current. Thus, the 
value of /2(0) is described by /2(C)) = (/%, — K n )fkim — 3/t n 4lm' where k p and K n represent the anomalous magnetic 
moments of the proton and neutron. In the case of neutron beta decay, we get /™^ p (0) = 1, 3™ ^ p (0) = F + D and 
f2~* p (0) = k p — k u , respectively. For the S° — > E + decay, the exact SU(3) symmetry predicts that /f^ s (0), fff^ s (0) 
and / 2 = '^ S (0) are identical to those of neutron beta decay. 

The experimental rate of the hyperon beta decays, B — > bW, is given by 



r = 



G 



60tt 



F ^(M B - M b f(l - 35)\V us \ 2 \f^ h (Q)[ 2 



sf^(0) 



/f-»(0) 



(6) 



where Gf denotes the Fermi constant. The ellipsis can be expressed in terms of a power series in the small quantity 
8 = (Mb — Mb) /(Mb + Mb), which is regarded as the size of flavor SU(3) breaking [l8|. The first linear term, which 
should be given by — 4(5[g2(0).gi(0)//i(0) 2 ]B^;,, is safely ignored as small as 0(8 2 ) since the nonzero value of the 
second-class form factor 512 should be induced at first order of the 8 expansion [13]. The absolute value of <7i(0)//i(0) 
can be determined by measured asymmetries such as electron- neutrino correlation 0, [l8| . Therefore a theoretical 
estimate of /i(0) is primarily required for the precise determination of |V^ S |. First of all, the value of /i(0) should 
be equal to the SU(3) Clebsch-Gordan coefficient up to the second order in 5/7(3) symmetry breaking, thanks to 
the Ademollo-Gatto theorem As the mass splitting among octet baryons is typically of the order of 10-15%, 

the expected size of the second-order corrections is a few percent level. However, either the size, or the sign of the 
second-order corrections are somewhat controversial among various theoretical studies at present as summarized in 
Table |H 

In the bag-model [2(| and quark-model calculations [2l], [52], flavor 5C/(3)-breaking effects on /i(0) are mainly 
accounted for wave- function mismatches between strange and non-strange quarks. Both models predict a small 
negative correction. On the other hand, the 1/N C expansion approach including SU(3) symmetry breaking up to the 
second order predicts a relatively large and positive correction. Recently, the full one-loop 0(p 4 ) calculation in heavy 
baryon chiral perturbation theory (HBChPT) was completed by Villadoro [24| . However, the author emphasized that 
the SU(3) version of HBChPT does not seem to be of help for the determination of /i(0). This is because a slow 
convergence of the chiral expansion is observed. It is also pointed out that a serious convergence problem is revealed 
by the inclusion of spin-3/2 decuplet degrees of freedom into the framework of HBChPT. Subsequently, the complete 
one-loop order result has been checked in a different regularization scheme, covariant baryon chiral perturbation theory 
(CBChPT) [25[. Both the size and the sign of the second order corrections are found to be different from results of 
HBChPT. The authors of Ref. [25[ have estimated partial corrections of 0(p 5 ) in HBChPT and then reconfirmed 
that the convergence behavior of SU(3) baryon chiral perturbation seems to be problematic as pointed out previously 
in Ref. [23| . Unlike the case of Ki 3 decays, the reliability of the chiral perturbation approach is questionable for the 
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hyperon decays. A model independent estimate of /i(0) is highly required to settle both the size, and the sign of the 
second order corrections on /i(0). 

The leading correction to the axial- vector coupling g\ (0) starts at first order in symmetry breaking, while flavor 
St/(3)-breaking effects to /i(0) are suppressed in first order by the Ademollo-Gatto theorem Therefore, sizable 
breaking corrections, which are of the order of 10% estimated from the mass splitting in the octet baryons, are to be 
expected in the ratio of <7i(0)//i(0). However, the current experimental precision is not enough to provide conclusive 
evidence of the violation of a two-parameter (F and D) fit based on the conventional Cabibbo model to ratios 
<7i(0)//i(0) measured in various hyperon decays @. As mentioned earlier, the 2° — » E + beta decay is highly sensitive 
to flavor SU(3) breaking since the ratio gi(0)//i(0) of this particular decay should be identical to that of neutron 
beta decay if the flavor S/7(3) symmetry is manifest. Therefore, flavor St/(3)-breaking effects may be easily exposed 
in the S° — » S + process. Indeed, the center-of-mass correction approach (26| and the 1/N C expansion approach [23| 
predict that the [<7i(0)//i(0)]e->e is smaller than the [fl , i(0)//i(0)]„_ p by 8-10% and 20-30% respectively. Such sizable 
breaking corrections could be distinguishable in experiment. However, the first and single experiment done by the 
KTeV collaboration at Fermilab showed no indication of flavor S[/(3)-breaking effects on <7i(0)//i(0), measuring 
fffi (0y A fOlla— >s — l-32±o;i7 [la]. The KTeV experiment reported no evidence for a nonzero second-class form factor 
g2 [161 ] within their experimental precision. The value of /2(0)//i(0) have been also measured in the KTeV experiment 
using the electron energy spectrum. Their observed value, [A(0)/A(0)]e->e = 3.8 ± 2.3, seems to be consistent with 
that of neutron beta decay as [A (0)/ A (0)1 

n —>p — K n — 3.706. Needless to say, its error is too large to discriminate 
either the exact ST/ (3) value or other theoretical predictions. See Table ITT1 where several theoretical predictions of 
the value [A(0)]e->£ are compiled. 

In this context, one would tend to conclude that the predictions of the exact ST/ (3) symmetry limit to hold better 
in the case of hyperons. Indeed, it is true that the mass splitting for hyperons is rather small compared to mesons. 
Nevertheless, as we will show from our lattice simulations, this is indeed not the case. 



TABLE I: Theoretical uncertainties of /i(0) for the H° — > E + transition process. 



type of result [/i(0)]eo_>e+ Reference 

bag model 0.97 [20] 

quark model 0.987 [21] 

quark model 0.976 [22] 

l/N c expansion 1.12±0.05 [23] 

full 0(p 4 ) HBChPT 1.009 a [24] 

full C(p 4 )+partial 0(p 5 ) HBChPT 1.004±0.026 [25] 

full Q{p 4 ) CBChPT 0.944+0.016 [25] 



The value is obtained by the iso-spin relation from that of the S — > S° transition process. 



TABLE II: Theoretical predictions of /2(C)) for the E° — ► E + transition process. For evaluations, we use current values of the 
anomalous magnetic moments of E and 3 baryons [27| . 



type of evaluation (Ref.) 


formula 




[/ 2 (0)] H o^ E+ 


[A(0)] s o^ B+ /[/ 2 (0)]„^ p 


exact SU(3) case 


Hp K.n 




3.706 


1 


Cabibbo model [9J 


M a +M s / n 
2M N Kn ) 




4.958 


1.338 


generalized CVC 






2.708 


0.731 


Sirlin's formula [2£] 


M-+M T 1 ,1 


)--|+^(« H0 + | Ks _) 


2.475 


0.668 


experimental value [16] 


N/A 




3.8 ±2.3 " 


1.03 ±0.62 



" A factor (Afg + Mj^)/Mw, equal to ~ 1.9048 is different from definitions of the /2 form factor adopted in Ref. |16( . 
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III. SIMULATION DETAILS 
A. Lattice set-up 

We have performed a quenched lattice calculation on a L 3 x T = 16 3 x 32 lattice with a renormalization group 
improved gauge action, DBW2 (doubly blocked Wilson in two-dimensional parameter space) gauge action [29|, [3(| 
at /3 = 6/g 2 = 0.87. The inverse of lattice spacing is about 1.3 GeV, set by the p-meson mass [3l|, yielding a 
physical volume of (2.4 fm) 3 . The spatial size 2.4 fm is large enough to accurately calculate the axial-vector coupling 
9 A = gi (0) [HI , which is one of the most sensitive observable to finite volume effects [3, IH, H3| • 

The previous quenched DWF studies by the RBC Collaboration reported that the residual chiral symmetry breaking 
of DWFs is significantly improved with a moderate size of the fifth-dimension. The residual quark mass for L s = 16 
is measured as small as m rcs ~5x 1CP 4 in lattice units (31|. Althou gh w e work with relatively coarse lattice spacing, 
a w 0.15 fm, good scaling behaviors of the light hadron spectrum [3l|. the kaon B-parameter Bk [H| and proton 
decay matrix elements [35[ are observed between at (3 = 0.87 (a « 0.15 fm) and 1.04 (a w 0.10 fm). Therefore, we may 
deduce that no large scaling violation is ensured for other observables as well in our DWF calculations. In Table IITT1 
some basic physics results are compiled from Ref. [3lj . 

In this study, DWF quark propagators were generated with three lighter quark masses m u( j = 0.04, 0.05 and 0.06 
for up and down quarks (60| and with two heavier quark masses m s =0.08 and 0.10 for the strange quark with 
L s = 16 and M5 = 1.8. We then take 5 different combinations between the up (down) quark and the strange quark as 
(m„ d ,TO s )=(0.04, 0.08), (0.05, 0.08), (0.06, 0.08), (0.04, 0.10) and (0.05, 0.10), which yield different S£/(3)-breaking 
patterns characterized by 5 = (Mb — Mb) /(Mb + Mb) for the B — > b process in the range of 0.009 to 0.028. Our 
results are analyzed on 377 configurations. Preliminary results were first reported in Ref. [36j (ol| . 

As mentioned earlier, the previous study of neutron beta decay with the same simulation parameters success fully 
yields a value of <7i(0)//i(0) as 1.212 ± 0.027, which just underestimates the experimental one by less than 5% [lj]. 
This success encourages us to study flavor SU (3)-breaking effects in the hyperon beta decays through a comparison 
between neutron beta decay and 2° beta decay. 

B. Mass spectra and dispersion relation 

In order to compute baryon masses or beta-decay matrix elements, we use the following baryon interpolating 
operator 

(Vx)ijk(t,P) =^e~ ip '*£abc[qa,i(yi,t)Cj 5 qb,j(y2,t)]q c , k (y 3 ,t) x 0(y x - x)0(y 2 - x)0(y 3 - x), (7) 

X 

where C is the charge conjugation matrix defined as C — jtjy and the index X S {B, b} distinguishes between the 
initial (B) and final (b) states. The superscript T denotes transpose and the indices abc and ijk label color and flavor, 
respectively. The superscript S of the interpolating operator 77 specifies the smearing for the quark propagators. In 
this study, we use two types of source: local source as </>(xi — x) = 5(xi — x) and Gaussian smeared source. Here we 
take xi = X2 = X3 = in our calculation. As for the Gaussian smeared source, we apply the gauge-invariant Gaussian 
smearing [37], HH with N = 30, u> = 4.35. Details of our choice of smearing parameters are described in Ref. [39l |. 

We construct two types of the two-point function for the baryon states. One interpolating operator at the source 
location is constructed from Gaussian smeared quark fields, while the other interpolating operator at the sink location 
is either constructed from local quark fields (denoted LG) or Gaussian smeared ones (denoted GG): 

C^ G (t-t sl . c ,p) - iTr{p+(?7|(t,p)^f(4rc,-p)} (8) 

with S = L or G. The projection operator V+ = Mp- can eliminate contributions from the opposite-parity state 
for |p| = [13, EH. It i s rather expensive to make the Gaussian smeared interpolating operator projected onto a 
specific finite momentum at the source location (t Brc ). However, it is sufficient to project only the sink operator onto 
the desired momentum by virtue of momentum conservation. Thus, the quark fields at the source location are not 
projected onto any specific momentum in this calculation. For the momentum at the sink location (£ s ink), we take all 
possible permutations of the three momentum p including both positive and negative directions in this study. 

All hadron masses are computed by using the LG-type correlators. We use the conventional interpolating operators, 
u^d (W75S) for the it (K) state, e a bc(Ua C^db)u c for the nucleon, £ a bc(UaCj5Sb)u c for the S state and £ a bc(s^ C^Ub)s c 
for the S state. All fitted values, which are obtained from the conventional single exponential fit for baryons (N, S 
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and 2) and the conventional cosh fit for mesons (w and K), are summarized in Tables IIVI and IVl Our simulated 
values of the pion mass range from 0.54 GeV to 0.67 GeV. 

The evaluation of the squared four- momentum transfer q 2 requires precise knowledge of the baryon energies Ex (p) 
{X = N, £, S) with finite momentum. This can be achieved by an estimation of the energy with the hel p o f the 
dispersion relation and the measured rest mass that can be most precisely measured. As we reported in Ref. [15| , the 
relativistic dispersion relation 



Ex(p) = ^Jp 2 +M x , (9) 

where p = (p x ,Py>Pz) with continuum-like momenta Pi — ^-^i (rii = 0, 1,2, ■ ■ -, (L — 1)), is indeed fairly consistent 
with the energies computed at least at the four lowest nonzero momenta: (1, 0, 0), (1, 1, 0), (1, 1, 1) and (2, 0, 0) in our 
simulations. It implies that our simulations do not much suffer from large 0(a 2 ) errors even at finite q 2 . We utilize 
such estimated energies instead of actually measured values in our whole analysis [CI- 



TABLE III: The residual mass m rcs , inverse lattice spacing (a^ 1 , set by the p meson mass), the renormalization factor of 
the axial-vector current (Za), the pion decay constant and the Kaon decay constant (Fk)- Those values are taken from 
Ref. [3l|. where simulations are performed on a 16 3 x 32 volume. 



Gauge action (/3) 


M 5 


L s 




a' 1 [GeV] 


ZAijrij = — m rcs ) 


Ft, [MeV] 


F K [MeV] 


DBW2 (0.87) 


1.8 


16 


5.69(26) xlO" 4 


1.31(4) 


0.77759(45) 


91.2(5.2) 


104.2(3.8) 



TABLE IV: Mass spectrum of non-strange hadrons (pion and nucleon) in lattice units. 



m ud M-k M N 



0.06 


0.5050(8) 


1.0821(42) 


0.05 


0.4617(9) 


1.0358(46) 


0.04 


0.4148(9) 


0.9869(50) 



TABLE V: Mass spectrum of strange hadrons (Kaon, E and H) in lattice units. 



m s m ud Mk My, M B 

0.08 0.06 0.5455(7) 1.1161(41) 1.1375(39) 
0.05 0.5257(8) 1.0895(43) 1.1210(39) 
0.04 0.5055(8) 1.0626(46) 1.1039(40) 
0.10 0.05 0.5652(8) 1.1201(46) 1.1741(40) 
(X04 0.5462(8) 1.0941(50) 1.1577(41) 



C. Three-point correlation functions 



We next define the finite-momentum three-point functions for the relevant components of either the local vector 
current (Sa) or t ne local axial current (3^ ) with the interpolating operators t\b and r\b for the B and b states: 

(rib(t',p')jZ(t,q)rj B (0,-p)) = g T a B ^ b { Pl p') X f(t,t',E B (p),E b (p')) + ■ ■ •, (10) 

where the initial (B) and final (b) states carry fixed momenta p and p' respectively and then the current operator has 
a three-dimensional momentum transfer q = p — p'. Here, Dirac indices have been suppressed. The ellipsis denotes 
excited state contributions which can be ignored in the case of if — t 3> 1 and t > 1. The ground state contribution 
of the three-point correlation function is described by two parts. The first part, Q^(p,p'), is defined as 

€' B ^ b (P,P') - (-i-y P ' + M b )O r a > B ^ h (q)(-i 7 -p + M B ), (11) 
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where 0^(q) corresponds to either Eq. @ or Eq. ([3]), and the factor f(t, t', Eb(p), Eb(p')) collects all the kinematical 
factors, normalization of states, and time dependence of the correlation function. The trace of Q^{p,p') with some 
appropriate projection operator V for specific combinations of T and a yields some linear combination of form factors 
in each T channel. On the other hand, all time dependences of the factor f(t, t' , Eb(p), Eb(p')) can be eliminated by 
the appropriate ratio of three- and two-point functions [42] 



B^b 



(*) = 



a 



V,B^b 



(t.p'.p) 



Cf, (^-sink ^src ) P ) 



Cfi (^sink ^)P)Cf, (t ^src; P')Cfo' < ^(^sink ^srcj P') 



Cu (tsmk ~~ t, p')C^ < ^(i — t SIC , p)C B (tsink — ^src; P) 



(12) 



where 



C^ b (t,p\p) = \TT{v(^(t siDk ,p')J^(t,ci)r^(t sl . c ,-p))} , 



(13) 



which are calculated by the sequential source method described in Ref. [141 ]. 

In this study, we consider the hyperon decay process B(p) — > 6(0) at the rest flame of the final (b) state (p' = 0), 
which leads to q = p. Therefore the squared four-momentum transfer is given by q 2 — 2Afb(^s(p)— -Mb) — (Afe — M,) 2 . 
The energies of the initial and final baryon state is simply abbreviated as Eb and E b , hereafter. In this kinematics, 
GaiPiP') is represented by a simple notation G^ip)- Then, the ratio (TT2")) gives the asymptotic form as a function of 
the current operator insertion time t, 



1 



1 



^/2M 2 Eb(Eb + M b ) 



(14) 
stc find 



in the limit when the Euclidean time separation between all operators is large, t s i n k > f > t src with fixed t 

^-sink- 

We choose particular combinations of the projection operator V and the current operator (r = V or A). Two 
types of the projection operator, V 1 — 7-+7t and V§ = 'P+757 2 are considered in this study. The latter operator 
implies that the z-direction is chosen as the polarized direction. We then obtain some linear combination of desired 
form factors from the following projected correlation functions: 



t n V.B^ b{q)} = Mb ( EB+MB 

t r y^ b{q)} = _ iqiMb 



E B 



M b + M, 



Mb fB - 

J2 



Tr{V§g^\q)} = - i£ijzqj M b [f 1 B ^ b (q 2 ) + ff^ b (q*)] 



V) 



Eb — Mb 



Mb fB ^ b 

~J3 



ft 



M h 



M b 



V) 



(15) 
(16) 
(17) 



for the vector currents J t and (i = x,y,z 



\Tr{V!g?' B ^ b (q)} 
\Tr{V£gf' B ^ b {q)} 



iq z M b 
M b 



Similarly, we get 
E B 



Mr + Mi 



Mb b- 

92 



Er 



(E B + M B )[g B ^ b (q 2 ) 



M F 



Mb + Mb 



M b B 

92 



M B + M b 



Mb B^b^ 



9?, 



V) 



M B + M b 



V) + 



(18) 



(19) 



for the axial-vector currents and (i — x,y,z). In this calculation, we use the four nonzero three-momentum 
transfer q = (n 2 = 1, 2, 3, 4). All possible permutations of the lattice momentum including both positive 
and negative directions are properly taken into account. All three-point correlation functions are calculated with a 
source-sink separation of 10 in lattice units, which is the same in the previous DWF calculations of the axial-vector 
coupling gA [H| and the weak matrix elements of the nucleon [l5| . 

Here, it is worth noting that the longitudinal momentum (q z ) dependence explicitly appears in Eq. (fP9|) due to our 
choice of the polarized direction. This fact provides two kinematical constraints on determination of the three-point 
functions in our calculation. First, there are two types of kinematics, q z ^ and q z = in the longitudinal component 
(i = z) of Eq. (|19j) . except for the case of n 2 = 3 where q z is always nonzero. Secondly, the transverse components 
(i = x or y) of Eq. (|19p are prevented from vanishing by the kinematics only if n 2 = 2 and 3, where two components 
of the momentum including the polarized direction (z-direction) are nonzero. 



8 



First of all, in Fig. [TJ we plot the dimensionless projected correlators of the vector part 

K v, B -, b j^g^M) (ok 
M b (E B +M B ) ' (2U) 

A v.B^ b = _I v lZ^!ggl!Ml fan 

s 3 . ^ i qi M b ' 1 J 

i— x,y ,z 

K v.B^ b = 1 / |Tr{Pggr- 5 ^ fc (g)} _ jTr{^g^ 5 - fc (g)} \ 

as a function of the current insertion time slice t for the S° — > S + process at (m u d,m s ) = (0.04,0.08) as typical 
examples. Good plateaus for all squared three-momentum transfer are observed in the middle region between the 
source and sink points. The quoted errors are estimated by a single elimination jackknife method. The lines plotted 
in each figure represent the average value (solid lines) and their 1 standard deviations (dashed lines) in the time-slice 
range 13 < t < 17. 

Similarly, we also define the dimensionless projected correlators of the axial-vector part as 

L M b (E B +M B ) ' 1 ' 

k A, B ^ = _i ( \^{nQ^ b m + k^{no^ b (q)} \ ; (24) 

T 2 ^ q z q x ' q z q y J ' 

A A,B^6 _ lTr{Vigf- B - b (q)} 

A ° - iq-M b ' (25) 

which are also evaluated from the plateau of the ratio (TT2")) . Fig. [5] is plotted for A^, which explicitly depends on the 
longitudinal momentum q z because of the chosen direction of the polarization. Two figures represent two types of 
kinematics, q z ^ and q z = 0. Solid and dashed lines are defined as in Fig. [JJ Good plateaus for all squared three 
momentum transfer q 2 are observed, similar to what we observe in the vector channel. 

The remaining two ratios and Aq are shown in Fig. [3J The upper figure is for A^, which is accessible only 
for q 2 = 2 and 3 in units of (2ir/L) 2 , where two components of the momentum including the polarized direction (z- 
direction) are allowed to be nonzero. Again, we observe reasonable good plateaus in the time-slice range 13 < t < 17. 
However, in the lower figure, the ratio A^ at two lower q 2 doesn't exhibit a clear plateau, while the reasonable plateau 
can be observed at two higher q 2 similar to other ratios. We count a short shoulder plateau in the time-slice range 
13 < t < 15 to take the average value of Aq at two smaller q 2 . It is worth mentioning that gi{q 2 ) and gz{q 2 ) are 
mainly determined by either A^ or A^, since contributions of A^ in Eqs. (|49|) and (|5"Tj) are numerically much smaller 
than others. The precise determination of the second-class form factor gi (q 2 ) may be affected by the poor plateau 
observed in A^ at two lower q. However, the subtraction procedure for the second-class form factors fz{q 2 ) and 
(?2(<Z 2 ), which will be described in the proceeding section (Sec. IVI TJ)) . may reduce the systematic error stemming from 
above mentioned issue. 



D. Renormalization 



In general, lattice operators receive finite renormalizations relative to their continuum counterparts since the exact 
symmetries of the continuum are usually realized only in the continuum limit a — > 0. Fortunately, the well-preserved 
chiral and flavor symmetries of DWFs [43|, [44|,|45[ make this task much easier than in the more conventional fermions. 
In this study, we use the vector and axial-vector local currents, which shares a common renormalization: Zy = Za, 
up to higher order discretization errors, 0(a 2 ) [46| . Therefore, we first focus on the vector renormalization. 

The vector form factors, especially in the precise determination of /i(0), require some independent estimation of 
Zy, the renormalization of the quark bilinear vector currents, 

[q f7a q f ,]™ = 4 S ' [977a9/'] latticc (26) 

where a subscript / denotes the flavor index. In this study, we need two vector renormalizations, Zy d and Zy s , for 
neutron beta decay and its SU(3) counterpart, the S° — > E + transition process. The former can be evaluated by the 
inverse of the forward limit of the n —> p vector matrix element because of lim q 2_ t0 {p\[u'jodY cn \n) = 1 in the present 
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FIG. 1: Relevant ratios of three- and two-point functions, Ao (top), Ag (middle) and Ay (bottom), for all possible three- 
momentum transfer q as a function of the current insertion time slice at (m u d,m s ) = (0.04,0.08). 
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FIG. 2: Relevant ratios of three- and two-point functions, h^(q z = 0) (top) and A^(g z 7^ 0) (bottom), for all possible 
three-momentum transfer q as a function of the current insertion time slice at (m u d, m s ) — (0.04, 0.08). 



calculation under the exact iso-spin symmetry (m u = mj). For the latter, this prescription is not directly applicable 
because of the presence of the flavor SU(3) breaking. However, we may calculate Z v s = Zy s (m u d, m s ) for m uc t ^ m s 
through the following relation: 

Z v s {m h m h ) = ^Z v s {m h mi)Z v s {m h ,m h ) (27) 

where mi and m/, (m; < m^) are simulated quark masses for either up (down) quark or strange quark. Zy s {mi,mi) 
and Zy s {mh,mi l ) correspond to the case of degenerate quark masses {m u d — m s ). Therefore, for an evaluation of 
those two vector renormalizations, we can utilize the relation, lim g 2_, (I] + |[u7os] Icn |S ) = 1, which is valid in the 
exact SU(3) symmetry limit. 

In the case of the flavor current (/ = /' in Eq. (|2"rj)) ). we had already observed that the relation Zy = Za is well 
satisfied in the chiral limit, up to higher order discretization errors 0(a 2 ) and neglecting explicit chiral symmetry 
breaking due to the moderate size of the fifth-dimensional extent L s . T his g ood chiral property of DWFs is known to 
be maintained even for the heavy- light vector and axial- vector currents [47| , which correspond to the extreme case of 
the flavor changing current (/ ^ /') in Eq. (|26[) . Therefore, in this study, we use the common renormalization given 
in Eq. (|2"T|) for both vector and axial-vector local currents. 
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FIG. 3: Relevant ratios of three- and two-point functions, Ay (top) and Aq (bottom), for all possible three-momentum transfer 
q as a function of the current insertion time slice at (m u d,m s ) — (0.04,0.08). 



IV. DETERMINATION OF /i(0) 
A. Scalar form factor fs(q 2 ) at q 2 = g, 



2 

max 



In the vector channel, only the time component of the vector current, namely the three-point correlation function 
jTr{V t Gt r ' B ^ t '(<?)} is prevented from vanishing at zero three momentum transfer |q| = by the kinematics This 



non-vanishing correlator gives the scalar form factor at specific four-momentum transfer as 



fs fc (9max) 



a V.B^b, 



0) 



(28) 



where g^ax = ~(Mb — M^) 2 [63|]. In Fig. [4j we plot the renormalized fsilmax) as a function of the current insertion 
time slice. Good plateaus are observed in the middle region between the source and sink points. The lines represent 
the average value (solid lines) and their 1 standard deviations (dashed lines) in the time-slice range 13 < t < 17. 
We stress that the statistical accuracy is less than about 0.5% even in the worst case {m u d, m s )=(0.04, 0.10). The 
obtained values of the renormalized fs(<irnax) as we ^ as ^ ne t> are one an d its renormalization factor Zy are summarized 
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7 2 ) 

1 max/ 



increases as m s deviates from m u( i, which was also observed 



in Table fVTl There is a tendency that the error of fs(q 2 
in the scalar form factor of K13 decays |3j • 

We should note that the renormalized /s(<7 max ) is exactly equal to unity in the flavor SU(3) symmetric limit, where 
/s(<7 max ) becomes /i(0). Thus, the deviation from unity in /s((7 max ) is attributed to three types of the SU(3) breaking 
effect: 1) the recoil correction (q 2 mx 7^ 0) stemming from the mass difference of B and b states, 2) the presence of 
the second-class form factor fs{q 2 ) and 3) the deviation from unity in the renormalized /i(0). Taking the limit of 
zero four-momentum transfer of fs(q 2 ) can separate the third effect from the others, since the scalar form factor at 
q 2 = 0, /s(0), is identical to /i(0). Indeed, to measure the third one is our main target. 
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FIG. 4: /s cn (g max ) as a function of the current insertion time slice. A source-sink location of three-point functions is set at 
[tsrc,i s ink] = [10,20]. The lines represent the average value (solid lines) and their 1 standard deviations (dashed lines) in the 
time-slice range 13 < t < 17. 



B. Interpolation to zero four-momentum squared 

The scalar form factor fs{q 2 ) at q 2 > is calculable with nonzero spatial momentum transfer (|q| 7= 0) We can 
make the q 2 interpolation of fs{q 2 ) to q 2 = together with the precisely measured value of fs(q 2 ) at q 2 = g^ ax < 0. 
First of all, we calculate the following combinations of two projected correlators ([20|) and ([2"T|) with nonzero three 
momentum transfer (|q| 7^ 0): 

/f-V) - EB - M » A^ b E B -M B VB ^ b 

Js yq ' Mb — Mi, M B -M b s y ' 

and then study the q 2 dependence of the scalar form factor. In Fig. [51 we plot the renormalized fs{q 2 ) as a function 
of four-momentum squared q 2 in physical units for {m uc i, m s ) = (0.04, 0.08) as a typical example. Either theoretically 
or phenomenologically, the q 2 dependence of fs(q 2 ) is not known due to the lack of knowledge of the second-class 
form factor f$. Our measured fs(q 2 ) up to at least q 2 < 1.0 GeV 2 exhibit a monotonic decrease with increasing q 2 . 
This observation is barely consistent with an expectation that fs{q 2 ) is dominated by fi{q 2 ), which is supposed to 
be the dipole form at low q 2 . 

In practice, the lack of the precise knowledge about the q 2 dependence of fs(q 2 ) is not a serious issue to determine 
/s(0) reliably. The simulated value of <? max is not far from q 2 = (65[. Therefore, /s(0) can be determined by a 
very short interpolation from g max , where we have very accurate data of /s(<7 max )- This allows us to expect that the 
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choice of the interpolation form does not affect the interpolated value /s(0) significantly. To demonstrate it, we test 
the monopole form 



fs(q 2 ) = 



fs(0) 

1 + a£V 



and the quadratic form 



/s(9 2 ) = /s(0)(l + A(V+A<V) 



(30) 



(31) 



for the q 2 interpolation of fs{q 2 ) to q 2 = 0. Indeed, both functional forms are adopted to evaluate /s(0) in a published 
work of the E~ — > n transition process [llj]. As shown in Fig. two determinations to evaluate fs(0) = /i(0) from 
measured points are indeed consistent with each other. All obtained values of /s(0) = /i(0) from both the monopole 
and quadratic form fits are summarized in Table IVHI Although the quadratic fit achieves the slightly smaller value 
of % 2 /dof ~ 0.04 than that of the monopole fit (~ 0.23), the highest q 2 point, which is not included in our fits, rather 
agrees with the monopole fit. We, therefore, do not use the results from the quadratic fit in the following discussion. 

Although it is hard to make a firm conclusion within the current statistical uncertainty, our measured values of 
/i(0) at the simulated points seem to receive small negative corrections of the SU(3) breaking. We then introduce 
the parameter of flavor SU (3) breaking, which is characterized by the measured mass difference between the 5 and £ 
states at the simulated points as 8 — (Ms — Ms) /{Ms + Ms). Our observed S[/(3)-breaking effect on /i(0), which 
corresponds to the deviation from unity, exhibits the quadratic dependence of this SU (3)-breaking parameter 6 in 
consistent with the Ademollo-Gatto theorem as shown in Fig. [6] Therefore, our results indicate that a sign of the 
second-order correction of the SU(3) breaking on /i(0) is likely negative. 




FIG. 5: Interpolation of fs(q 2 ) to q 2 — for {jn uc i,m s ) = (0.04,0.08). Open circles are fs(q 2 ) at the simulated q 2 . The solid 
(dashed) curve is the fitting result by using the monopole (quadratic) interpolation form, while the open diamond (square) 
represents the interpolated value to q 2 — 0. 



C. Chiral extrapolation of /i(0) 



In order to estimate /i(0) at the physical point, we perform the chiral extrapolation of /i(0). As described earlier, 
/i(0) can be parameterized as 



A(0) = 1 + A/, 

where A/ represents all SU(3) breaking corrections on /i(0). We then introduce the following ratio: 

A/ 



R M {M K ,M„) = 



{Mi-M*r 



(32) 
(33) 
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FIG. 6: The deviation from unity on /i(0) increases as the 5(7(3) breaking <5 increases. All measured values of /i(0) exhibit 
quadratic scaling with respect to the S'C/(3)-breaking parameter, which is suggested by the Ademollo-Gatto theorem. The 
dashed (quadratic) curve is a guide to the eye. The physical point is represented by a filled square. 



where the leading symmetry-breaking correction, which is predicted by the Ademollo-Gatto theorem, is explicitly 
factorized out As reported in Ref. [3], the remaining dependence of the SU(3) breaking is hardly observed 

within the statistical errors. To examine the quark mass dependence, we plot -Ra/, which is evaluated by the monopole 
form for the q 2 interpolation, as a function of M\ + M 2 in Fig. [71 There is no appreciable dependence of simulated 
quark masses within the statistical errors. This observation is consistent with what we observe in FigJHl where all 
measured values /i(0) at different simulated quark masses exhibit quadratic scaling with respect to the SU (3)-breaking 
parameter S. We then consider a linear fit in terms of M|- + M 2 : 

Raj (Mr, M„) =A +Ai- {M 2 k + M 2 ). (34) 

A dashed line in Fig. [7] corresponds to the chiral extrapolation using the linear form (|3"4")l . We obtain the extrapolated 
value of Ra / at the physical point as 

i?A/(M£ hys , Mf ys ) = -0.22(24) in (GeV) -4 (35) 

by employing result from the monopole form for the q 2 interpolation. We finally quote 

[A rcn (0)] M = 0.987(19) (36) 

at the physical point. By combining with a single estimate of |Kxs/f^ S (0)| = 0.216(33) from the KTeV experi- 
ment [la ], we obtain 

\V US \ = 0.219(27) cxp (4) thoory , (37) 
which is consistent with the value obtained from Ki 3 decays and the CKM unitary predicted value [l], 0, [1] • 



V. DETERMINATION OF 3i(0)//i(0) 

A. Ratio gi(q 2 )/f.s{q 2 ) at q 2 = ql ax 

Let us consider the ratio of A^' B ~ >h and Aq' B ~ >b , which are both accessible at zero three-momentum transfer 
q| = in the axial-vector and vector channels respectively. From this ratio at |q| = 0, we can evaluate the value of 
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FIG. 7: Chiral extrapolation of -Ra/- Open circles (diamonds) are calculated with m a =0.08 (0.10). The extrapolated point at 
the physical point is represented by a filled square. 



TABLE VI: Results for Zy tt (m„,[,m s ), /s"(?Lx) and /s°(9„„), which are evaluated in the region of [t mill ,t ma x] = [13, 17]. 





m ud 




/if** (<?max) 


/^"('/max) 


0.08 


0.06 


0.8144(7) 


1.2286(24) 


1.0010(15) 




0.05 


0.8128(8) 


1.2310(37) 


1.0010(27) 




0.04 


0.8112(10) 


1.2327(58) 


1.0005(44) 


0.10 


0.05 


0.8168(7) 


1.2240(46) 


1.0003(35) 




0.04 


0.8152(10) 


1.2248(65) 


0.9990(51) 



TABLE VII: Results for [/s en (0)] H ^ E = [/i° n (0)] s „, E , which are evaluated by the q 2 extrapolation with the monopole and 
quadratic functional form. 





m ud 


[/r n (o)] s ^ E 

monopole 


quadratic 


0.08 


0.06 


1.0004(16) 


1.0006(16) 




0.05 


0.9986(27) 


1.0001(28) 




0.04 


0.9985(44) 


0.9989(44) 


0.10 


0.05 


0.9965(35) 


0.9972(35) 




0.04 


0.9938(50) 


0.9947(51) 




phys. point 


0.9868(191) 


0.9897(192) 



9l(q 2 )/fs(q 2 ) atq 2 



Q 2 

imax 3 



max) 



A,B^b 



V,B^b 




(q = Q) 

(q = 0) 



(38) 



where Q^ax = —(Mb — M\,) 2 . We note that this ratio is exactly equal to <7i(0)//i(0) in the flavor SU(3) symmetric 
limit, which corresponds to that of n — ► p. In Fig. [5J we plot the ratio as a function of the current insertion time 
slice. Good plateau behaviors of <h(9max)//.s(9max) are observed in the middle of region between the source and sink 
points. Results of ih(<Zmax)//s(9max)' which are summarized in Table IVlnl are averaged over appropriate time slice 
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range 14 < t < 16. 

In our exploratory study (36| . we proposed the following double ratio: 

Af- s (q^o).A ^(q=o) / gr s (g, 2 naj-^ s (gLJ \ / f gr p (on r , Q1 

A y ^ s (q = 0) • A^(q = 0) U^tea*) - Sf?**^))/ l/r(0)J ' ^ J 

where 5 = . Since this double ratio is exactly equal to unity in the flavor SU (3) symmetric limit, the deviation 

form unity exposes flavor SU (3)-breaking effects in the S — > E decay. As shown in Fig. [5J the double ratio exhibits 
good plateau behavior slightly above unity in the time slice range 13 < t < 17. The deviation from unity becomes 
large as increasing the size of the flavor SU (3)-breaking, which is characterized by the size of S as listed in Table IVTTTI 
The observed deviation indeed contains three types of the SU(3) breaking effect, similar to what we explained 
for fsiQmax)- Here, we note that <7i(0)//i(0) receives the first-order corrections since the axial-vector form factors 
are not protected by the Ademollo-Gatto theorem. Therefore, we expect that the flavor SU (3)-breaking observed in 
the double ratio could be dominated by the leading symmetry-breaking correction on [<?i(0)//i (0)]e->s- The reasons 
as follows: 1) The q 2 dependence of form factors at q 2 — qf nliK , which is proportional to 6 2 , can be involved in the 
second-order corrections as /f~* S (i?max) = /f^ S (0) + 0{8 2 ) and gf^ S (<?max) = 5i^ S (0) + 0(S 2 ). 2) The nonzero 
value of the second-class form factors and g 2 starts from the first-order corrections. These contributions in g% and 
fs are involved in the second-order corrections as well. As a result, the double ratio is expressed by 

[Si(0)//i(0)] n _ >p 

where the first term is responsible for the leading first-order correction. The deviation from unity observed in the 
double ratio may be able to exhibit the size of the leading SU (3)-breaking correction on g\ (0) / /i (0) for small 5 (6(| . 
As listed in Table I Villi the observed size of flavor SU (3) breaking effects is indeed comparable of the size of the 
leading order 0(6) [361 ] and glows linearly with increasing the parameter S. 

We simply perform the linear fit in two mass combinations M\ + M 2 and M|- — M 2 on Rd as 

R D = B + B 1 ■ (M 2 K + Ml) + B 2 • (Ml - M%). (41) 

We then obtain Rd = 1.022(31) at the physical point. A sign of the observed corrections to [<7i(0)//i(0)]e->x: 
seems to be opposite to model predictions from the center-of-mass correction approach [26| and the X/N c expansion 
approach [23J. However we recall that the observed corrections less than a few percents are too small to justify 
neglect of the second-order corrections in our analysis since the natural size of the flavor SU(3) breaking is around 
10% [67|]. To make a firm conclusion, we have to removes systematic uncertainties induced by neglecting both the 
recoil corrections and the presence of the second-class form factors. 



TABLE VIII: Results for ffi(?max)//s(9niax)) which are evaluated in the region of [t m in,imax] = [14, 16]. The double ration Rd 
and the SU(3) breaking parameter S are also listed. 





m ud 


3l(<?m a x)//s(gmax) 


Rd 


x _ M s -M E 


0.08 


0.06 


1.221(13) 


1.0042(28) 


0.0095(3) 




0.05 


1.219(15) 


1.0061(55) 


0.0142(5) 




0.04 


1.218(17) 


1.0099(103) 


0.0190(7) 


0.10 


0.05 


1.225(13) 


1.0114(75) 


0.0236(7) 




0.04 


1.223(15) 


1.0151(126) 


0.0283(10) 


phys. point 




1.0217(309) 


0.04996 



B. Results for <?i(0)//i(0) and ffi(0)//i(0) 

The recoil corrections are removed by taking the limit of considered form factors to q 2 = 0. In similar to the case of 
fs(q 2 ), the value of <h(0)//s(0) can be evaluated by the q 2 interpolation of gi{q 2 ) / fs{q 2 ) to q 2 = 0. The form factor 
9i^ b (q 2 ) at <? 2 > can be calculated by the projected correlator (f23|) with zero longitudinal momentum (q z = 0) but 
nonzero transverse momentum (q x or q y ^ 0) [68j: 



9 ^\q 2 )=k A L B -\q z =Q). 



(42) 
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FIG. 8: 9i(gmax)//s° n (9max) as a function of the current insertion time slice. The lines represent the average value (solid lines) 
and their 1 standard deviations (dashed lines) in the time-slice range 14 < t < 16. 
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FIG. 9: The double ratio Rd as a function of the current insertion time slice. A source-sink location of three-point functions is 
set at [t S rc,isink] = [10, 20]. The lines represent the average value (solid lines) and their one standard deviations (dashed lines) 
in the time-slice range 13 < t < 17. 



One can calculate the ratio of gi(q 2 )/ fs(q 2 ) with Eqs. ([28]) and (j42|) . In Fig. [TU1 we plot the ratio of gi{q 2 )/ fs(q 2 ) 
as a function of q 2 . We consider two types of the interpolation form, the monopole and quadratic forms, the same 



18 



TABLE IX: Results for <7i (0)//i (0) , which are evaluated by the q 2 extrapolation with the monopole and quadratic functional 
form. 





m ud 


5i(0)//i(0) 
monopole 


quadratic 


0.08 


0.06 


1.220(13) 


1.220(13) 




0.05 


1.217(15) 


1.217(15) 




0.04 


1.215(17) 


1.215(17) 


0.10 


0.05 


1.222(13) 


1.222(13) 




0.04 


1.220(15) 


1.220(15) 




phys. point 


1.205(27) 


1.206(27) 



as in the case of fs(q 2 ). The highest q 2 data point is not included in our fits. All obtained values of <?i(0)//s(0) = 
<7i(0)//i(0) from both the monopole and quadratic form fits are summarized in Table ITXl Again we observe that two 
determinations to evaluate <?i(0)//s(0) from measured points are fairly consistent with each other. Therefore, we do 
not use the results from the quadratic fit in the following discussion, the same in the case of /s(0). 

Next, we perform a linear fit in two mass combinations JW| + M 2 and M|- — M 2 for the values of <7i(0)//i(0) and 
then obtain [<7i(0)//i(0)]e-»e = 1-205(27) at the physical point. It may be compared with its St/(3)-syrnmetric value 
of 1.191(49), which corresponds to the chiral extrapolated value of [gi(0)//i(0)] n _ >p to the physical point by using 
the simple linear fitting form in terms of M 2 . These results give [<7i(0)//i(0)]e^e = 1.016(31) x [ffi(0)//i(0)] w -» p . 

After the q 2 dependence is taken into account, the deviation from unity is now reduced by 0.6% from the double 
ratio Rd- If the conventional assumption 52(0) = is adopted here similar to usual experimental analyses, the flavor 
St/(3)-breaking found in <7i(0)//i(0) tends to be tiny. Although it does not conflict with the Cabibbo theory, the 
following alternative interpretation is still not ruled out. It is possible that a relatively small first-order correction on 
<7i(0) is accidently canceled out in the form g\ = g\ — 8gi by an opposite contribution stemming from the large and 
positive value of (72(0) such as <?2(0)/.gi(0) ~ 1. As we will discuss in next section, this is indeed the case. Therefore, 
we have to subtract the contribution of the second-class form factor properly in order to estimate the true size of the 
first-order correction on <7i(0)//i(0). 

A complete analysis requires information of the g% form factor. The individual form factors in Eqs. @ and ([3]) can 
be determined at finite |q|. Then, the value of the ratio of 52 (9 2 ) / 9i(q 2 ) at zero momentum transfer are obtained 
by an appropriate q 2 extrapolation. We finally obtain <72(0)/<7i(0) = 0.677(177) at the physical point. See the next 
section for details. 

After the subtraction of the 52(0) contribution, we obtain [ffi(Q)//i(0)]s-»s = 1.248(29) at the physical point. It 
implies that the relatively large and positive value of <?2(0)/<7i(0) has induced a shift of the value of g\ (0)/fi (0) toward 
the exact St/(3)-symmetric value. Finally, we obtain 



ffi(O) 
/i(0) 



= 1.051(35) x 

♦E 



3i (0) 



/i(0) 



(43) 

n — >p 



at the physical point. The deviation from unity is increased by 3.5% in comparison with the unsubtracted case. 
Although the size of the St/(3)-breaking corrections on <7i(0)//i(0) is still relatively smaller than the expected size 
evaluated from the mass splitting among the octet baryons (~ 10%), the similar size of the flavor St/ (3) breaking in 
the S + — ► n decay was also reported in Ref. [131 ] . 

In Fig. 111! we summarize our result and the experimental values combined with predictions from the center-of- 
mass correction approach [2(| and the X/N c expansion approach 23]. Although the experimental data is not yet 
sufficiently precise to determine either the size, or the sign, of the SU (3)-breaking corrections, our result suggests that 
the symmetry-breaking correction is likely small but positive. It is worth mentioning that the sign of our observed 
corrections is opposite to the model predictions, but in agreement with that of the E + — -> n decay measured in 
quenched lattice QCD [l]|. 
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FIG. 10: Interpolation of gi{q 2 )/ fs(q 2 ) to q 2 = for (m ud ,m s ) = (0.04,0.08). Open circles are gi{q 2 )/ fs(q 2 ) at the simulated 
q 2 . The solid (dashed) curve is the fitting result by using the monopole (quadratic) interpolation form, while the open diamond 
(square) represents the interpolated value to q 2 — 0. 
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FIG. 11: Comparison among model predictions, experimental data (KTeV-FNAL) and our lattice result for the ratio of 
[ffi(0)//i(0)]H-E and its SU(3) counterpart [ffi(0)//i(0)]„^ p . 



VI. RESULTS FOR OTHER FORM FACTORS 



A. Computational Method 

The kinematics of |q| 2 = allows only a particular combination of the projection operator (V) and the Lorentz 
index of the currents (a) in either vector or axial- vector channels fl4| . However, in the case if spatial momentum 
transfer q is nonzero, all three-point correlation functions defined in Sec. IIII CI are calculable. Therefore, three form 
factors at finite |q| can be obtained individually by solving simultaneous linear equations. For the vector channel, the 
simultaneous linear equations are given by 

= . (44) 



Eb-Mb E B -M b 



M B +M b M B +M b 

E B -M b Eb + Mb 

M B +M h M B + M b 

1 
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One gets each form factor by inverting the above equations [69| as 
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(45) 
(46) 
(47) 



Similarly, three-point correlation functions of the axial- vector part are also described by the following simultaneous 
linear equations, 
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and then each individual form factor is given by inverting the above equation as 



(48) 



9^" 



9i-» 



(9 2 ) 



M B + M b 

2Mb 
Mb + Mb 

2Mb 
Mb + Mb 

2Mb 



A 



A.B^b 



A L (ft 



a A,B—tb/ 



0) 
0) 



Mb - 
M B + M b 

A,B^b 



A 



A.B^h 



B 



A 



o 



-A 



M b 

A,B^b 



Mf, ^A,B^b 



0) + A 



A.B^b 



Mb 

Eb + Mb 



A 



A,B^b 



(49) 
(50) 
(51) 



Here we remark that as described previously, we utilize all possible permutations of the lattice momentum including 
both positive and negative directions and adopt four nonzero values of three- momentum transfer q = (n 2 = 1, 
2, 3, 4). Here, it should be reminded that the ^-direction is chosen as the polarized direction in this study. This 
fact makes the analysis of the axial-vector channel more complex than the vector channel. Indeed, the longitudinal 
momentum (q z ) dependence explicitly enters in A^. Accordingly, there are two types of kinematics, q z ^ and 
q z = in the three-momentum transfer, except for the case of n 2 = 3 where q z is always nonzero. In other words, 
A r ' B_> (q z = 0) at n 2 = 3 can not be calculated directly. However, instead, we may evaluate A L ' B ^ b (q z = 0) at 



n 



3 from the longitudinal correlator A^' B ^ b and the transverse correlator A^,' B ^ b by using a relation 



A 



A,B^b 



(ft = 0) = A 



A,B- 
L 



'(ft o) 



M b {E B + M B ) 



-A 



A,B^b 



(52) 



which is easily read off from Eq. (| 19|) . Note that Af / (q z ^ 0) are always calculable at finite |q| 



B. Second-class form factors: f'i(q 2 ) and g-2{q 2 ) 
1. Subtraction method 



Figure [T2l shows the ratios of fz{q 2 ) / fi{q 2 ) (left panel) and 92{<l 2 ) / 9i{<1 2 ) (right panel) as a function of four- 
momentum squared q 2 for (m u d,m s ) — (0.04,0.08). Open circles represent the values measured for 5 — > £ at the 
simulated q 2 . Although we observe non- negligible values of the second-class form factors, it is still questionable 
whether nonzero signals correspond to a pure effect from the flavor SU (3) breaking. Indeed, the same analysis for 
the case of the exact ST/ (3) limit (n — > p), yields comparable values of the second-class form factors to those of 
S — > E. Lower and upper triangle symbols correspond to the cases of n — > p with m u d — 0.04 and 0.08. The lighter 
m u d calculations yield central values closer to those of S — > £ with larger statistical uncertainties, while results from 
the heavier m u d calculations are marginally consistent with zero values for both fa and 52 form factors within 1-2 
standard deviation. Although it seems that observed nonzero values of the second-class form factors suffer much from 
large statistical fluctuations, we are also concerned about another possibility. 

For the case of the f% form factor, a nonvanishing contribution even in the exact Sf/(3) limit stems from a subtle 
difference of Ag and A^ correlators, which is possibly due to the discretization error [70]. Although there is no 
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such clear correspondence in the case of the gi form factor, the discretization error may equally cause a systematic 
uncertainty on the determination of the gi form factors as well. 

To control both statistical fluctuations and systematic uncertainties, we subtract the measured values of the second- 
class form factors in the SU (3) limit calculation from those of 3 — > E as 

[/rV)] H ^ - m 2 )}-^ - [/3(9 2 )v , (53) 



[fl5 u V)] B ^ E = [92(q 2 )] 



3^S 



2Mjv 
M s + Af s 
2M N 



92(q 



(54) 



where a factor (Me + Me)/(2Mjv) is accounted for the mass difference between the 3 — > E and n — > p decays, which 
is determined by simulated masses [7l| . 

Although the q 2 value for the 3 — > E transition differs from that of the n — > p transitions at the same three- 
momentum transfer, above subtraction is simply performed at every given three-momentum q. There are choices to 
set a reference value of the subtraction, since the second-class form factors for the n — ► p transition are supposed to 
vanish with any value of m u d- However, as mentioned above, the lighter m uc t calculations provide larger statistical 
uncertainties on the second-class form factors than those of 5 — > E. In Eqs. (|53[) and ([54]), we adopt the single 
reference values of [fz{q 2 )} n ^ p and [g2(q 2 )]n^ P evaluated at m u d = 0.06. Our chosen value of m u d corresponds to the 
heaviest m u d mass in all combinations of (m tt d,m s ) for the 3 —> E calculation in this study. 

Results from the subtraction method are also plotted in Fig. [T3] as filled circle symbols. It is clearly observed that 
statistical errors are significantly reduced after such subtraction due to a strong correlation between those two form 
factors, while center values are slightly shifted to zero. Non-vanishing signals of both and g2 form factors turn out 
to be more pronounced. The subtraction method can expose the real £[/(3)-breaking effect with better statistical 
precision. 
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FIG. 12: The subtraction method can expose non-zero values of the second-class form factors, which corresponds to the 
Sf/(3)-breaking effect, with better statistical precision. 



2. Extrapolation to zero four-momentum squared 

We next examine the q 2 dependence of the ratios fi uh (q 2 )/ fi(q 2 ) and g2 Uh (l 2 ) / 9i(q 2 ) ■ In Fig-HSl we show the case 
of (m u d,m s ) = (0.04,0.08) as typical examples. One can easily observe that both ratios do not yield any strong q 2 
dependence at least in our simulated range of 0.25 to 0.93 GeV 2 . It is worth mentioning that there is no theoretical 
prediction of the q 2 dependence of the second-class form factors. Therefore, we simply adopt the linear or quadratic 
extrapolation with respect to q 2 to estimate /3(0)//i(0) and <?2(0)/<7i(0). All evaluated values of /3(0)//i(0) and 
#2 (0)/<?i (0) from two functional forms are summarized in Table IXl The extrapolated values from both determinations 
agree well with each other within their errors. As a result, the extrapolated values are not significantly affected by 
the specific fitting form adopted for describing the q 2 dependence of form factors. 
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FIG. 13: The ratios }3{q 2 ) / }i{q 2 ) (left panel) and g2(q 2 ) / gi(q 2 ) (right panel) as functions of q 2 as in the case of (m u d,m B ) = 
(0.04,0.08). The solid and dashed lines represent linear and quadratic fits for three lowest q 2 data, respectively. 



3. Chiral extrapolation 

In this subsection we show the results of /3(C)) //i(0) and 52(0) /<?i(0) obtained by two types of the chiral extrapo- 
lation. In the first method, which we call method-A, we take into account the fact that the second-class form factors 
fs and <72 vanish in the exact SU (3) symmetry limit and nonzero values are induced by the first-order corrections. 
We then introduce the following ratios for both /3(0)//i(0) and <72(0)/<7i(0): 



d (M M s /3(0)//i(0) 



R 92/ gi {M K ,M V ) 



g2(0)/.9i(0) 



(55) 
(56) 



where the leading first-order corrections of the flavor SU(3) breaking are explicitly factorized out. In Fig. [21 we 
show the quark mass dependences of Rfe/h and R g2 / gi as functions of M|- + M%. The ratio Rf 3 /f 1 reveals the mild 
dependence of Mjj + M 2 , while rather strong dependence of M\ + M 2 appears in the ratio R 92 / gi - In either cases, 
obtained data is well described by the linear fitting form [72| : 



R h/h (Mk,M v ) or R 92/gi (M K , M„) = A a + A x ■ {M 2 K + M*). 



(57) 



Here, the values of /a(0) and 52(0) given by the simplest linear fit in q 2 are used for calculating /3(C)) //i(0) and 
52(0)/si(0). We then get the extrapolated values of Rf 3 /f 1 and R g2 / gi at the physical meson masses as 



Rh/h(Mf ys ,M^ s ) = 0.61(42) in (GeV)- 2 , 
R 9 ,/gAM l K yS ,Mt yS ) = 3.02(89) in (GeV)- 2 , 



(58) 
(59) 



which finally provide the values / 3 (0)//i(0) = 0.137(94) and ff 2 (0)/si(0) = 0.677(177). 

In an alternative method indicated by method-B, we may perform a linear fit in two mass combinations + M 2 
and M\ — M% 



/s(0)//i(0) or fla(0)/ 5 i(0) = Bo + B 1 ■ {M 2 K + M 2 ) + B 2 ■ (M 2 K - M 2 ) 



(60) 



We then obtain / 3 (0)//i(0) = 0.147(60) and g 2 (0)/gi(0) = 0.450(110). Both method-A and B provide consistent 
results. Although the errors in the latter approach are relatively smaller than that of the former, the former leads to 
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a much smaller value of % 2 /dof than the latter. Therefore, we quote the values obtained from the method-A for our 
final values at the physical point: 



/a(0) 



/i(0) 
92 (0) 



9i (0) 



3^S 



= 0.137 ±0.094, 



= 0.677 ±0.177, 



(61) 
(62) 



which show firm evidence for nonzero second-class form factors in the 3 — ► E + beta decay. It should be reminded that 
the KTeV experiment reported no evidence for a nonzero second-class form factor gi [161 ]. measuring g2(0)//i(0) = 
-0.89 ± 1.05, which corresponds to 3 2 (0)/si(0) ~ -0.73 ± 0.89 [z|[- 
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FIG. 14: Chiral extrapolation of Rf 3 /f ± (left panel) and R g2 / gi (right panel). Symbols are defined as in Fig. [6] 



TABLE X: Results for ]/3(0)//i(0)] H ^ E and [32(0)/(/i(0)] H ^ E , both of which are evaluated by the q 2 extrapolation with the 
linear and quadratic functional form. 







]/ 3 (0)//i(0)] s ^ E 




faa(0)/ffi(0)]a^ s 




m s 


m ud 


linear 


quadratic 


linear 


quadratic 


0.08 


0.06 


0.028(13) 


0.011(22) 


-0.017(26) 


0.027(52) 




0.05 


0.045(13) 


0.037(22) 


0.048(26) 


0.069(49) 




0.04 


0.064(20) 


0.067(36) 


0.131(38) 


0.124(76) 


0.10 


0.05 


0.062(22) 


0.040(35) 


0.030(43) 


0.088(83) 




0.04 


0.078(26) 


0.066(42) 


0.107(49) 


0.143(92) 


(From method-A) 


phys. point 


0.137(94) 


0.241(173) 


0.677(177) 


0.414(372) 


(From method-B) 


phys. point 


0.147(60) 


0.193(110) 


0.450(114) 


0.344(232) 



C. Weak magnetism form factor f2{q 2 ) 



In Fig. [T5l we show the weak magnetism form factors /;? (<? 2 ) for 3 — > S (left panel) and n — > p (right panel) 
at (m u d,m s ) = (0.04,0.08) as functions of four-momentum squared q 2 . The form factors plotted here are scaled by 
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the renormalization factor Zy determined in Sec lIVl To determine /2 Cn (0), we adopt two functional forms for the q 2 
dependence of f^il 2 )'- the dipole form 



/TV) 



/ 2 ren (o) 



(l + X^q 2 ) 2 '' 



and the quadratic form 



/rV) = /r(o)(i + a<V + a< 2 Y)- 



(63) 



(64) 



The former functional form is phenomenologically adopted in the nucleon's electromagnetic form factors, which are 
related to the weak nucleon form factors under the CVC hypothesis on the weak processes. 

Figure [TBI shows that both functional forms are equally fitted to data of /| on (g 2 ) in either case of 3 — > E or n — > p 
decay processes. However, there appears to be a sensitivity of the choice of the fitting form in extrapolated values at 
q 2 = 0. This is simply because our simulated q 2 points are not close enough to q 2 = 0. In Ref. [l5j], we have studied 
the q 2 dependence of the weak nucleon form factors at low q 2 down to about 0.1 GeV 2 , which is accessible with the 
spatial extent L = 24. Indeed, in our previous study, the weak magnetism form factor for n — > p is observed to be well 
described by the dipole form [15]. To make this point clear, in the right panel of Fig. [TBI we also include results from the 
larger volume [L = 24) for a comparison. The steep raising behavior of L = 24 data as q 2 decreases clearly favors the 
dipole form. We, therefore, do not use the results from the quadratic fit in the following discussion. All extrapolated 
values with both determinations are summarized in Table IXIl As for the n —> p transition, the extrapolated values of 
/2(0) by using the dipole fit are summarized in Table IXIT1 together with other relevant quantities. 

We also perform the chiral extrapolation of /""(O) through a linear fit in two mass combinations M|- + M 2 and 
M 2 , the same as Eg. ([601) . We obtain [/""Wis— »£ = 3.30 ± 0.24 by employing result from the dipole form for 



the q 2 extrapolation. Finally, we compare the ratio [/2(0)//i(0)]e 
then observe an order of 16% breaking effect as 



>s to its SU(3) counterpart [/2(0)//i(0)]„_ >p and 



/a(0) 



/i(0) 



1.16(11) 



/i(0) 



(65) 



n — >p 



which implies the violation of the exact SU(3) relation. It is worth mentioning that [/ 2 (0) / /i (0)]s_ > j: tends to be bigger 
than [/2(0)//i(0)]n-»p, while the /2(0)//i(0) values evaluated from both the generalized CVC hypothesis [12J and 
Sirlin's formula [28| yield opposite results as previously shown in TableQj] In addition, the observed size of the deviation 
from unity in Eq. (frJS")) may also be compared with the Cabibbo-model prediction as (M= + M^)/(2M N ) = 1.338 0], 
which corresponds to a factor accounted for the mass difference between the S — > E and n — > p decays. We again 
observe a 15% deviation from the Cabibbo model. 



TABLE XI: Results for renormalized [/2(C))] s _ s is evaluated by the q 2 extrapolation with the dipole and quadratic functional 
form. 





m u d 


[/2 Cn (0)] 2 ^ E 

dipole 


quadratic 


0.08 


0.06 


3.64(11) 


3.25(12) 




0.05 


3.60(13) 


3.20(14) 




0.04 


3.54(15) 


3.14(16) 


0.10 


0.05 


3.63(11) 


3.23(12) 




0.04 


3.57(13) 


3.17(14) 




phys. point 


3.30(24) 


2.92(25) 



D. Induced pseudo-scalar form factor gz (q 2 



We next show the q 2 dependence of the ratio of the induced pseudo-scalar form factor 173 (q 2 ) and the axial- vector 
form factor g\{q 2 ) at (m u( i,m s ) = (0.04,0.08) in Fig. [TBI It is phenomenologically known that the q 2 dependence 
of flsfg 2 ) in the low q 2 region are well fitted by the pion-pole dominance (PPD) form, gy,{q 2 ) — 4M 2 / g 1 (q 2 )/(q 2 + 
M 2 ) [H, [5(| HH- A s shown in Ref. |15| . the PPD form provides a good description of the q 2 dependence of the 
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FIG. 15: The renormalized weak-magnetism form factor f2(q 2 ) as functions of four- momentum squared q 2 for H° — > E + (left 
panel) and n — > p (right panel) at (m u d,m s ) = (0.04, 0.08). For the q 2 extrapolation, we apply two types of fitting form. The 
solid (dashed curve) is the fitting result by using the dipole (quadratic) form, while the open diamond (square) represents the 
extrapolated value to q 2 = 0. The highest q 2 data point is not included in our fits. In the right panel, results from the larger 
volume (L = 24) [l5|] are also included as open triangles for a comparison. 



TABLE XII: Summary for the n — > p transition process. The values of /2(0)//i(0) determined by the dipole form for the q 2 
extrapolation, while the values <?3(0)/<;i(0) are evaluated by the pion-pole fit. 





ffi(0)//i(0) 


/a(0)/A(0) 


ffs(0)/ 5 i(0) 


0.10 


1.240(12) 


3.80(8) 


12.06(32) 


0.08 


1.226(11) 


3.71(10) 


13.52(38) 


0.06 


1.214(14) 


3.57(13) 


15.35(49) 


0.05 


1.209(18) 


3.47(17) 


16.73(59) 


0.04 


1.202(24) 


3.35(24) 


18.68(79) 


phys. point 


1.191(49) 


3.18(29) 


151.0(6.9) 



nucleoli's induced pseudo-scalar form factor measured in quenched DWF simulations at least at low q 2 . On the other 
hand, there is no direct experimental information in the case of gs{q 2 ) for any other hyperon decays. 

Under the partial conserved axial- vector current (PCAC) hypothesis [z3], an extension of PPD to AS* = 1 decays 
predicts that the induced pseudo-scalar form factor for the hyperon beta decays, such as the S — > X decay, the K 
meson pole, instead of the pion pole. In the generalized PPD form, the ratio of gz(q 2 )/ gi {q 2 ) for the S — > £ transition 
is given by a simple monopolc form 



53 (q 2 



LSi(9 2 ) 



PPD 



(Afe + Afs) 2 
q 2 + M 2 K 



(66) 



where a monopole mass corresponds to the Kaon mass. 

In Fig. [TBI we include the predicted q 2 dependence (dotted curve) evaluated by the generalized PPD form with 
simulated baryon masses {Mm, Ms, M~) and meson masses (M w , Mk)- Three larger q 2 data points are quite close 
to predicted curves in either 3 — > X or n — > p, while the lowest q 2 data points are underestimated in comparison with 
the PPD prediction. To extrapolate the value of g^(q 2 ) / ' gi{q 2 ) to zero four-momentum transfer, we first adopt the 
monopolc form, which is inspired by the PPD form. In Fig. 1161 the dashed curves are fitting results by the monopole 
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form with two parameters A3 and A 3 ' , 



where corresponds to a monopole mass. The obtained values of <?3(0)/<h(0) are significantly smaller than 

those of the PPD prediction as (M s + M s ) 2 /M|- for the S -> E transition and 4M^/M^ for the n — ► p transition. 
The resulting monopole masses are also bigger than the expected Kaon and pion masses. This substantial deviation 
from the PPD form is mainly caused by an unexpected reduction of the lowest q 2 data points. 

In Ref. [Hj], we have reported that the lowest q 2 point in the nucleon's induced pseudo-scalar form factor at L = 16, 
which is now utilized in our current calculations, may suffer from the finite volume effect. Indeed, the q 2 dependence 
of the results obtained from the larger lattice (L = 24) was well fitted by the pion-pole structure. In this context, 
we may have an alternative way to evaluate the value of 33(0)/gi(0) by a one-parameter fit in the monopole form 
with the monopole mass constrained by simulated Mk or . We simply refer to such a fit as the "Kaon-pole fit" or 
the "pion-pole fit" respectively. All extrapolated values of gz{q 2 ) / ' gi{q 2 ) to zero four-momentum transfer with both 
determinations of monopole and Kaon-pole fits are listed in Table [XTTTl 

The solid curves in Fig. [TC] are given by the Kaon- and pion-pole fits. For a justification of this analysis, see the 
right panel of Fig. [TH] In this figure, four data points (open triangles) obtained from the larger lattice (L = 24) are 
additionally included. One can easily see that these data points quite follow the solid curve, which is determined by 
the pion-pole fit of the lowest three q 2 data points obtained from the lattice size of L = 16. It is found that resulting 
values of <?3(0)/gi(0) are still slightly smaller than the PPD values in either case of S — > £ or n — > p. The similar 
quenching was observed in our previous detailed study of the nucleon's induced pseudo-scalar form factor using the 
larger lattice (L — 24) [15]. It is worth mentioning that the size of this quenching for S — > E is similar to that 
of n — > p. Therefore, ratios of [53(0) / <7i(0)]e->£ and [53(0) / <?i(0)] n _v P exhibit remarkable consistency with the PPD 
prediction as shown in Table IXIVI Even if we adopt the monopole form for the q 2 extrapolation, the resulting ratios 
still barely agree with the PPD values. Our results strongly suggest that the following relation is well fulfilled at least 
in the simulated region 

" ff 3 (0) 
.5i(0) 

As for the chiral extrapolation of 53 (0) jg\ (0) , the linear fit in two mass combinations M% + M 2 and Mjj — M 2 (the 
same as method-B in Sec. I VI B~3|) was utilized in Ref. 13;]. However, this extrapolation doesn't take into account the 
expected large quark-mass dependence of g$ (0) j g\ (0) in the vicinity of the chiral limit like a divergent 1 /M 2 term for 
n — > p or a 1/Mjj term for 5 — > E. This implies that the extrapolated values should be considerably underestimated 
especially for the case of n — > p. For the ratio of [<?3(0)/<?i(0)]e->£ and [03(O)/</i(O)] n _>p at the physical point, the naive 
chiral extrapolation indeed yields a large value of 0.67(16) for the monopole fit or 0.52(4) for the Kaon(pion)-polc fit, 
which should be compared with the PPD value of M\{Ms + M s ) 2 /(4M| r M|) = 0.1430 at the physical point. This 
result is clearly contradicted with the finding expressed by Eq. (|68[) fulfilled in the simulated region. 

The simple linear fit in two mass combinations M% + M 2 and M\ — M 2 is instead applied to the ratio of <?3 (0) / g\ (0) 
and its PPD value, which has a very mild quark-mass dependence in either case of S — > E or n — > p. The value 
of [s , 3(0)/pi(Q)]h— >s at the physical point is evaluated from the extrapolated value of this ratio with a multiplicative 
factor of the physical PPD value. We obtain [5 3 (0)/ffi(0)] H ^ s = 21.58(98) and [g 3 (0) / gi(0)] n ^ p = 151.0(6.9) at the 
physical point for the Kaon(pion)-pole fit. The ratio of those values, which are determined to be 0.1429(2), correctly 
reproduces the PPD value. This determination is rather phenomenological. However, it is hard to perform the chiral 
extrapolation of <?3(0)/<7i(0) without any assumption within our limited data sets. Thus, instead of quoting any final 
value, we would like to stress that the expected relation between [<73(0)/<7i(0)]e_»e and [<73(0)/<?i(0)] n _,.p as in Eq. ([68]) 
is confirmed in our simulations. 



VII. SUMMARY 

In this paper, we have studied flavor ST/(3)-breaking effects in the hyperon semileptonic decay, S° — > Y> + lvi using 
quenched DWF simulations with the lattice size L 3 xT = 16 3 x 32. The spatial extent La w 2.4 fm was large enough 
to calculate all six form factors describing the beta-decay matrix element without a serious finite volume effect on 
the axial-vector coupling <?i(0). From phenomenological point of view, the significance of this subject is twofold: (1) 
to extract the element V us of the Cabibbo-Kobayashi-Maskawa mixing matrix from the AS = 1 decay process, and 



g 3 (g 2 ) 

5i (q 2 ) 



1 + A 3 V 



(67) 



a— >y: 



2M 



N 



Ml 
Ml 



93(0) 
5i(0). 



(68) 



i> — ► (1 
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FIG. 16: The ratio gs{q 2 )/ gi{q 2 ) as functions of four-momentum squared q 2 for E° — > S + (left panel) and n — > p (right panel) 
at (m u d, m s ) = (0.04, 0.08). Symbols are defined as in Fig. 1151 The highest q 2 data point is not included in our fits. 



TABLE XIII: Ratio [ff3(0)/<7i(0)] n _, B is evaluated by the q 2 extrapolation with the monopole and Kaon-pole form. 









fos(0)/$i(0)]a^s 






m u d 




monopole 


Kaon pole 


0.08 


0.06 




12.08(87) 


14.40(42) 




0.05 




12.36(98) 


14.90(46) 




0.04 




12.63(1.15) 


15.44(52) 


0.10 


0.05 




11.77(86) 


13.95(42) 




0.04 




11.97(1.00) 


14.34(47) 


phys. point 




16.73(2.51) 


21.58(98) 


TABLE XIV: Comparison of measured ratios of [ff3(0)/g , i(0)] = _ lS 


and [ffa(0)/5i(0)] n ^ p to their PPD value 


given by M%(M B + 


M S ) 2 /(4M 2 N M 2 ). 










m a rn u d 


From monopole fit 


From Kaon(pion)-pole fit 


PPD value 


0.08 0.06 


0.957(17) 


0.938(6) 


0.930(1) 


0.05 


0.925(33) 


0.891(10) 


0.878(2) 


0.04 


0.882(58) 


0.827(15) 


0.811(3) 


0.10 0.05 


0.881(43) 


0.834(14) 


0.819(4) 


0.04 


0.836(68) 


0.768(18) 


0.751(5) 



(2) to provide vital information to analysis of the strange quark fraction of the proton spin with the polarized deep 
inelastic scattering data. Our particular choice of the S° — > S + decay process is highly sensitive to the flavor SU(3) 
breaking, since this decay process is nothing but the direct analogue of neutron beta decay under the exchange of the 
down quark with the strange quark. The SU (3)-breaking pattern observed in this study is summarized in Table IXVI 
The vector form factor at zero four- momentum transfer, /i(0) is protected by the Ademollo-Gatto theorem against 
corrections at first order in symmetry breaking. However, a sign of the second-order correction is somewhat con- 
troversial among various theoretical studies at present. Our estimate of renormalized [/i(0)]e— >s at the physical 
point from quenched lattice QCD simulation is 0.989(19), which indicates that the second-order correction on /i(0) 
is likely negative. This leads to the closer value of \V US \ to the value obtained from K13 decays. Although both the 
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TABLE XV: Summary of the SC/ (3)-breaking pattern observed in the H° — > E + decay. Here, A§u(3) is defined as ([X]b->£ — 
[X]n->p)/[X]n_> p for each quantity X. We also evaluate the scaled values of A§u(3) by the expected size of the leading-order 
corrections. 



X 


A-SU(3) 


A su(s) /(2S) n 


n 


/i(0) 


-0.013(19) 


-1.3(1.9) 


2 


MQ) 


+0.16(11) 


+1.6(1.1) 


1 


Mo) /Mo) 


+0.137(94) a 


+1.4(9) 


1 


Si(0)//i(0) 


+0.051(35) 


+0.51(35) 


1 


fla(0)/ffi(0) 


+0.677(177) a 


+6.8(1.8) 


1 



Because /3(C)) = and 92(0) = for n p, [/a(0)//i (0)]3_>j; and [g2(0)/gi(0)]g_>£ are instead quoted respectively. 



1/N C expansion analysis and the full one-loop 0(p 4 ) calculation in HBChPT favor positive corrections, our observed 
tendency for the SU(3) breaking correction agrees with predictions of quark models and CBChPT up to complete 
0{p-). 

The leading correction of the flavor SU(3) breaking to g\ (0) starts at first order in symmetry breaking. Although 
sizable breaking corrections, which is the order of 10% estimated from the mass splitting in the octet baryons, 
is expected, we found relatively small and positive correction to <?i (0)//i (0) as [g\ (0)//i (0)]e->e = 1.051(33) x 
[51(0) / /i(0)] n _»p in contrast to the model predictions where large and negative correction is preferable. Unfortunately, 
the first and single experiment of the S° — » S + decay done by the KTeV Collaboration is not yet sufficiently precise 
to determine either the size, or the sign, of the SU(3) breaking correction to <?i(0)//i(0). 

The advantages of lattice QCD studies of the hyperon beta decay are further demonstrated in determinations 
of the other beta-decay form factors, while it is difficult to determine each form factor separately in experiments. 
Especially, information of the second-class form factors gi is required since linear combinations of g\ (0) and 52 (0) are 
actually measured in the experiments from the Dalitz plot that reflects the electron- neutrino angular correlation [l8| . 
Furthermore, the nonzero value of the weak electricity form factor gi as well as that of the induced scalar form 
factor /3 is the direct evidence of the SU(3) breaking effect in the hyperon beta decays. We obtain the ratios of 
[ff2(0)/5i(0)] S ->E = 0.68(18) and [/ 3 (0)//i(0)]h-e = 0.14(9). Although both values are roughly comparable to the 
expected size of the leading first-order corrections of the flavor SU (3) breaking, the former is much larger than the 
latter. 

A remarkable observation is that a relatively small first-order correction to <?i(0) is accidentally compensated for 
flavor S'[/(3)-brcaking effects on <?i(0) with an opposite contribution due to the relatively large and positive value of 
172(0). This may suggest why the conventional analysis of the hyperon beta decays based on the Cabibbo hypothesis 
works well, though the effects due to the SU(3) breaking observed in the octet baryon masses and magnetic moments 
are expected to considerably affect the axial-vector part of the weak matrix elements. 

We have also found that the weak magnetism /2(C)) receives positive corrections of order 16 % for the flavor ST/ (3) 
breaking, measuring [/2(0)//i(0)]e->£ = 1.16(10) x [/2(0)//i(0)] n _ p . Our result is not in agreement with either the 
generalized CVC hypothesis or the Cabibbo-model prediction. On the other hand, as for the induced pseudo-scalar 
form factor 33, the generalized PPD form, which is extended even in AS — 1 decays under the strong assumption 
of PCAC, provides a good prediction of the ratio of [03(O)/0i(O)]e-»e and [g 3 (0)/gi(0)], wp at the physical point as 
Mj,(M s + M s ) 2 /(4:M^MZ) = 0.1430. This indicates that the large SC/^-breaking effects on g 3 (0)/#i(0) is attributed 
to the Kaon(pion)-pole structure of the 173 form factor. 

In this study, we have focused only on the specific beta-decay process, S° — > E + , However, if the flavor SU(3)- 
breaking pattern observed here would commonly appear in other beta-decay processes, our results call for non- 
negligible SU (3)-breaking effects in all hyperon beta decays. We also believe that the quenched approximation is not 
problematic for the determination of flavor SU (3)-breaking effects in the hyperon beta decays in similar to what was 
observed in calculations of K13 decays [H, Nevertheless, the simulation with dynamical 2+1 flavor quarks is an 

important future direction to be explored for full knowledge of the SU (3)-breaking pattern in the hyperon beta decays. 
Especially, in order to settle the signs of the leading order correction on /i(0), /2(0) and <7i(0)//i(0), more extensive 
lattice study is required. We plan to extend the present calculation to include other relevant hyperon beta-decay 
processes such as £~ — * n and A — > p using dynamical Nf = 2 + 1 flavor DWF lattice configurations generated by 
the RBC and UKQCD Collaborations [H, [5J| . Such planning is now underway [55| 
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Appendix A: Scalar function fs 



Let us consider the matrix element of the divergence of the vector current: 



(b(p')\d a V a (0)\B(p)) = u b (p')[i(S>-ll)f^\q 2 ) 
= (M b -M B )f^ b (q 2 )- 



M 



Mb + Mb 

" 2 A^V) 



f^ b (q')]u B (p) 



M h 



u b (p')u B (p)- 



(69) 



Here, we have used the Dirac equation for both initial (B) and final (b) baryon states, (ip + Mb)ub{p) — Ub(p')(ifl 
Mb) = to get from the first line to the second line. Combined with Eq. (fj|, one finds the following relation 



(b(p')\d a V a (0)\B(p)) = (M b - M B )fE^ b {q 2 )ub{p')u B {p)i 



(70) 



where an overall factor Mb — Mb on the right hand side is responsible for the current conservation when the flavor 
SU(3) symmetry is exact (Mb = Mb)- 



Appendix B: Other parametrization of the baryon weak matrix element 

Instead of the standard parametrization of Eqs. ^ and ([3]), the following equivalent form [56| is more useful to 
derive all of Eqs. (fP^ - (fl~9|) . which are considered at the rest flame of the final (b) state (p' = 0): 

/nV ( \ rB~*br 2\ . ■ f2^ b (l 2 ) . ■ fz^iQ 2 ) n 
a (?) = 7a A (q)+W a — — TT+^aTT ~TT '> ( 71 ) 



Mb + Mb Mb + Mb ' 

' M B + M b iqalb M B + M b ' 



o a (q) = 7a755i (g) + *Pa75T7 — nrr + ^7577 — rrr- ( 72 ) 



The two sets of form factors are connected by 



A 3 ^ V) = A 3 - V) + A 5 - V), g? ^ V) = g?- b (f) fe^^V), 

fj- b (q 2 ) = 2fi^(q 2 ), ~g^ b (q 2 ) = 2g^ b (q*), (73) 

./^ V) = /J^V) A^ V), ~9^ b {l 2 ) = 9^ b {q 2 ) ~ .9^ V), 
One can easily check above relations using the Gordon identity. 
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